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Abstract 

We derive boundary conditions at interfaces (contact discontinuities) 
for a class of Lagrangian models describing, in particular, bubbly flows. 
We use these conditions to study Kelvin-Helmholtz' instability which de- 
velops in the flow of two superposed layers of a pure incompressible fluid 
and a fluid containing gas bubbles, co-flowing with different velocities. We 
show that the presence of bubbles in one layer stabilizes the flow in some 
intervals of wave lengths. 



1 Introduction 

Many mathematical models of fluid mechanics are derived through the approx- 
imation of the solution of boundary value problems for Euler equations. For 
example, equations for bubbly flows are derived as an approximation of the 
solution of a complex free-boundary problem describing the motion of the mix- 
ture of water and gas bubbles (lordanski (1960), Kogarko (1961), Wijngaarden 
(1968)). Green-Naghdi's model describing wave motions of a liquid layer of fi- 
nite depth in a shallow water approximation with account of dispersion eflFects 
is another example (Green et al (1974), Green & Naghdi (1976)). Due to the 
fact that averaging procedures and asymptotic expansions have been used in 
the derivation, it is not obvious to decide what boundary conditions are natural 
for these systems of equations. In this paper, a Hamiltonian formulation of the 
problem in (t, x)-space is proposed which allows one to find boundary conditions 
for a general class of models. This class includes a model of bubbly fluid and 
dispersive shallow water. We use the boundary conditions to study the stability 
of two co-flowing layers of bubbly and pure fluids. In absence of gravity and 
capillarity, it is well known that Kelvin-Helmholtz' instability develops for any 
wave lengths in the flow of two superposed layers of pure incompressible fluids 
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(sec, for example, Drazin & Rcid (1981)). We prove that the presence of bub- 
bles in one layer can produce a stabilizing effect: the flow becomes stable with 
respect to perturbations of some wave lengths. 



2 Governing equations and conditions on mov- 
ing interfaces 

2.1 Variation of Hamilton's action 

Here we calculate the variation of Hamilton's action for a special class of La- 
grangians. Usually this procedure is used only for derivation of governing equa- 
tions. We obtain below not only governing equations but also boundary condi- 
tions at inner contact surfaces. 

Let us consider the Lagrangian of the form: 

L = L(J,-,z) (2.1) 

where z = ^ ^ ^ = (z*) , (i = 0, 1, 2, 3), t is the time, x is the space variable, 

J = ^ ^ . p is the fluid density and u is the velocity fleld. Let us calculate 

the variation of Hamilton's action in the case when the 4-D momentum J verifies 
the equation of continuity: 

Div J = ^ + div{pu) = (2.2) 
Hamilton's action is defined by 

Ldz (2.3) 



Jq 



where is a material domain. Considering a smooth one-parameter family of 
virtual motions 

z = *(Z,£), *(Z,0)=<p(Z) 

(Z stands for the Lagrangian coordinates, e is a small parameter at the vicinity 
of zero and z — ipCZ) is the real motion), we define the virtual displacements 
C(Z) and the Lagrangian variations 53{Z) by the formulae: 



C(Z)=^*(^'^) 



de 



5m - ^'^^'^^ 



e=0 ^£ 



(2.4) 



£=0 



Due to the fact that Z = <p~^(z), we consider the variations as functions of 
Eulerian coordinates and use the same notations ^(z) and (5J(z) in Eulerian 
variables. Hamilton's principle assumes that ^(z) = on the boundary dVl of 

n. 
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-Div (A^SJ) - Div (A^) S3 



dL 



and finally 

5a — I tr 
Jn 



.dL 
dJ 



d3\ dC 
dz J dz 



dc 



dL 



L - — j) Div C + Div {A^SJ) 



dL 



dz 



tr 



dJ 



dz 



dC 
dz 



dL 



dL 



L - — J + (DivA^) J] Div C + Div {A^6J) + C 



dz 



Let us denote 



Then 



dL 



dJ 



dz 



5a ^ [ — - Div [ J - A' — + IL - J) L 



dz 



(2.8) 



+ / f J - A^^ + (L - J) h] C + N^A'^SJ 
Jon \ dz J 

with N'^ — {—Dn,ri'^), Dn denotes the surface velocity and n is the space 
unit normal vector. Virtual displacements vanish at the boundary dfl and the 
surface integral is zero. The volume integral yields the equations of motion in 
conservative form as in Gavrilyuk & Gouin (1999). 

In the case when fluid tensorial quantities are discontinuous at the inner interface 
E, expression ()2.8|) becomes 



Sa 



dL 



Div J - 



5J 

d^ 



L-K' J) h 



(2.9) 



NM J - A 



.d3_ 
dz 



L-K' J) h ] C + '^^A^SJ 



where the jump through E is denoted by [ ]. 
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2.2 Governing equations and boundary conditions 



We explicit the governing equations and the inner boundary conditions for the 
Lagrangian 



L = ^p\u\'^ ~W (^p,pj , where () 



dt dt 



u^V() (2.10) 



Such a Lagrangian appears in the study of wave propagation in both shallow 
water flows with dispersion and bubbly flows (a complete discussion of these 
models is given in Gavrilyuk & Teshukov (2001)). We get 



d3_ 



I dp^ dp_ \^ 

dt dx 

\ dt dx / 



dL 
9J 



|ur dW IdW , ,T T IdW , n 
— (Vp) u, u --— (Vp) 



2 dp p dp 



P dp 



Matrix (12.71) becomes 



dL 



5z 



dW 
dp 

dW 
V dp ' 



O3 



We get 



Div{A^ 



u 



d_ dW 
dt\dp, 

6W 1 dW 



, dW \ n 

div u , 0^ 

\dp ) 

1 dW 



-^-^ + -— (Vp) u, - -— (Vp) 



P dp 



J = p (n^ u - Dn) , N'^A^ 



P dp 

dW , rp „ N 7 

u-Dn) , 0^ 

dp 



(2.11) 



We study the case when S is a contact surface and, consequently, J = 0, 
N^A^ _ qT -yyg -|.-|^^-|. -j-jj^g surface integral in (|2.9p vanishes if the boundary 
conditions : 

[L-K^ J] = [p] =0 (2.12) 

are fulfilled with p = L — J. Vanishing the volume integral in (|2.9p and 
using relations (j2.1ip . we obtain the governing equations in the form 



dp 
dt 



+ div{pu) = 0, 



(2.13) 



^ + div (puu^+pl) = 0, 
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A complete set of boundary conditions at the contact interfaces for the class of 
models considered is: 

n^u-Dn = 0, [p] = (2.14) 

Equations (|2.13p have been obtained earlier in Gavrilyuk & Shugrin (1996) by 
means of the classical method of Lagrange multipliers. Notice that this method 
did not give jump conditions (I2.12p . System (I2.13P is reminiscent of Euler 
equations for compressible fluids. The term p in equations of motion (I2.13P 
and boundary conditions (j2.14|) stands for the pressure. Nevertheless, p is not a 
function of density as in the case of barotropic fluids: it depends also on material 
time derivatives of the density. Thus, p is not usual thermodynamic pressure. 
In general, it implies that the boundary conditions are not just consequences 
of conservation laws. For example, for Korteweg-type fluids where the stress 
tensor depends on the density gradient, the boundary conditions contain normal 
derivatives of the density at interfaces (Seppecher (1989), Gouin & Gavrilyuk 
(1999)). As we have shown here, the boundary conditions at inner interfaces in 
bubbly fluids do not depend on normal derivatives of the density: they involve 
only tangential derivatives. 

If the flow domain is shared between two domains consisting of pure ideal fluid 
and bubbly fluid, all previous calculations are still valid and boundary conditions 
(|2.14p can be extended to such a case. 

2.3 Exact statement of the problem 

Now, we propose to use conditions (|2.14p to study Kelvin-Helmholtz' instability 
of a parallel flow between two rigid walls of a layer of bubbly fluid in contact 
with a layer of pure incompressible fluid. 

Consider the 2-D flow of two superposed finite layers of inviscid fluids in the 
channel: — oo < x < +oo, < y < Hq. The first one is governed by the Euler 
equations : 

div(ui) = 0, — i + (uTv) ui H = 0, -oo <x<oo, <y < hit, x) 

ot Pw 

. (^-^^^ 

Here pio is the constant density of pure fluid, Ui — (ui, vif^ is the velocity field 
and pi is the pressure. 

In the domain h{t, x) < y < Ho governing equations (|2.13p for the second fluid 
are 

^ + divip2U2) = 0, (2.16) 
ot 

— + (U^ V) U2 + ^ 0, 

OT P2 

P2^P2{^~^A — ]~dlv[ — U2 ]]-W, W = W{p2,P2) 




dt 



'2 
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Here p2, U2 = {u2,V2)'^ and p2 are the average density, velocity and mixture 
pressure, respectively; W = W{p2,P2) is a given potential. For the flow of 
compressible bubbles having the same radius, potential W has the form (see, 
for example, Gavrilyuk (1994) and Gavrilyuk & Teshukov (2001)): 

W{p2,P2) = P2 (cgSgipg) ^ 2mipiR^R^ (2.17) 

where Eg is the internal energy of the gas in bubbles, R is the bubble radius, pg 
is the gas density, pi — const is the density of carrying phase, Cg = const is the 
bubble mass concentration, n is the number of bubbles per unit mass. Potential 
W is the difference between the internal energy of gas and the kinetic energy 
of fluid due to radial bubble oscillations. The bubble radius and the bubble 
density are functions of the average density p2 '■ 

^rf.i(l-i^), ,,..,(1-1^)- 

3 n \p2 pi J \p2 pi J 

It can be shown that governing equations (|2.16p with potential (|2.17p coincide 
exactly with classical equations of bubbly fluids (Kogarko (1961) and van Wi- 
jngaarden (1968)). In particular, the role of "equation of state" 

8W 

P2^P2-, W, W = W{p2,P2) (2.18) 

0P2 

plays the Rayleigh-Lamb equation which governs the radial oscillations of a 
spherical bubble: 

RR+^R ^I{pg-P2), Pg^pl^ (2.19) 
I Pi ^ dpg 

The equivalence between p.l9p and (I2.18P with W given by ()2.17|) has been 
proved, for example, in Gavrilyuk (1994). Dissipation-free model p.l6|) - p.l7p 
assumes that the sliding between components is negligible. Moreover, this model 
is valid only when the volume fraction of bubbles is very small. Notice that a 
detailed description of the bubble interaction has been recently done by Russo 
& Smereka (1996) and Herrero, Lucquin-Desreux & Perthame (1999) in the case 
of rigid bubbles, and Smereka (2002), Teshukov & Gavrilyuk (2002) in the case 
of compressible bubbles, by using a kinetic approach. 

At the rigid walls y — Q and y = Hq, the vertical components of the velocity 
are equal to zero 

vi {t, 0) ^V2{t,x,Ho) =0. 

In accordance with the previous Section, we prescribe the following boundary 
conditions at the contact interface y — h(t,x): 

ht+Uiha:=Vi, ht+U2hx^V2, Pi ^ P2 (2.20) 
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3 Linear stability problem 
3.1 Linearization 

Consider the following main parallel flow of two fluids: 

ui ~ uio — const, vi — 0, Pi = Pq — const, < y < /iq, = const 

(3.1) 

U2 = U2Q = const, V2 = 0, P2= Pa = const, p2 = P20 = const, ho < y < Hq 

(3.2) 

Small perturbations denoted by ' : 

wi=wio + Ui, vi^Vj^, pi=pq+Pi, 

P2=P2Q+P2, U2=U20+U2, V2 = V^, P2 ^ PO + P2 

satisfy the linearized system of equations 

uix + viy=0, pio-DiUi = 0, pio-Dit;! +J3i,y = 0, < y < ho (3.3) 

D2P2+P20{u2x+'^2y) = 0> P20D2U2+IP2X = 0> P2oD2v' +P2y =0, ho < y < Ho- 

(3.4) 

We use the notations 

Di = -k: +uoi — , 1 = 1,2 
at ax 

In (EH) 

P2 = a^P2 + ^^^2^2 

2 X ,2 5^^/ 

" = P20^-^(P20,0), = -P20 2"('°20,0). 

Here a is the equilibrium sound velocity, and 6 is a characteristic wave length 
depending on the bubble radius and gas volume fraction. We suppose that 

— -(P20,0) = 

dp2ap2 

This condition is obviously fulfilled for potential (|2.17p . The coefficients a and 
b calculated in equilibrium R = Ro and p2 = po are : 

„2 _ dpg 1 , 2 _ 1 Pi (or. 

- dr^^'>NoP2o' - A^NoRoP20 ^ ' 

For dilute mixtures we can simplify these expressions by replacing p2o by pi : 

2 dpg , , 1 2 1 ^0 

a = — 3~(^o)-r7 — , b = - 



with 



dr No pi AttNoRo 3ao 
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Here Nq = p2on is the number of bubbles per unit volume, R — Rq is the 

4 

equilibrium bubble radius, uq = tqNq is the volume fraction of gas, tq = —ttRq 

3 

is the bubble volume and pg = pg (r) is the gas pressure in bubbles expressed as a 
function of bubble volume. The equilibrium sound speed a is usually small with 
respect to the gas sound velocity. Expressions (|3.5p can be obtained directly 
after linearization of Rayleigh-Lamb's equation for bubbles (I2.19p . 
Boundary conditions (|2.20p a,t y — ho are 

/ / / / / / / / 

hf + uioh,j. = v-^^, hf + U2oh^—V2, Pi — P2 (3-6) 

Let us consider the normal modes of the linear problem (|3.3p . (|3.4p . (|3.6p : 

Ul — Ui{y) eyip{i{kx — Ljt)), — ikVi{y) eyip{i{kx — cut)), (3.7) 

Pi — Pi{y) exp(i{kx ~ Lut)), h ^ H exp(i{kx ~ cut)) 
P2 = -^2(2;) exp(i(A;x - ujt)), = U2{y) exp(i(fca; - ut)), 
V2 = ikV2{y) exp{i{kx — uit)), P2 = Piiy) exp{i{kx — uit j) 

By substituting into the linearized system we get the system of equations for 
unknown amplitudes 

Ui + Viy = 0, pio(uio-c)C/i+Pi -0, (3.8) 

-pi0fc^(Ul0 - C)Vl + Ply = 0, 

{U20 - c)i?2 + ^20(^^2 + V2y) = 0, P20("20 - C)U2 + P2 = 0, 
-P20k\u20 - C)V2 + P2y =0, P2 - - 6'fc'(w20 ^ cf) i?2 

with the following boundary conditions 

yi(0) = 0, y2(i?o)=0, Pi{Ho) = P2{Hol Yl^ = Y^ = H 

uio - c M20 - c 

Here c = — is the phase velocity. 
k 



From equations (|3.8p we obtain the eigenvalue problem for pressure amplitudes: 
Piyy -k^Pi^O, 0<y< ho, (3.9) 
P2yy ~ X^k^P2 =0, ho<y< Ho, 

Plyjho) ^ P2y[ho) 
Plo(wiO-c)2 P20[U20~CY 

where 



Piiho) = P2iho), „ \.2 - „ .,2 ' Piy{0) = P2y{Ho)^Q 



,2 -, {U20-Cf 



d?- — \P'k'^(u2o — c)^ 
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It follows from (|3.9p that the eigenvalues c are solutions of the equation 
thikho) A th{Xk {Hq - ha)) 



pio{uiQ - c)^ P2o(w20-c)2 



==0 (3.10) 



Next we assume that both layers are thin: k {Hq — ho) <C 1, kho <C 1, and 
dispersion relation (jS.lOp is simplified into 

(Ho I 1 ^ , h, ^ ^3^^^^ 



P20 \("20-c)2 - 62fc2(M20 - c)2/ pio(uio - c) 
3.2 Study of dispersion relation 

In case M20 = uio, equation (|3.1ip reduces to the quadratic equation 



(1 + A)a^ - ((1 + A) b^k^ + 1) (m20 - = 0, with A 



{Ho - ho)pio 



which has only real roots. It means that the flows with equal velocities are 
stable. 

Let us consider the general case when U20 ^ uiq. Equation (|3.1ip can be rewrit- 
ten as a polynomial of fourth degree. The stability needs that all four roots are 
real. To study this problem, it is convenient to rewrite the dispersion relation 
(|3.1ip in the form 

F{z)+A = (3.12) 

with 

F{z) = {l-Nz)^(^l--^y (3.13) 

U20 — c ok a a 

Here M is similar to the Mach number, d is the dimensionless length of pertur- 
bation wave. Obviously, z is real if and only if c is real. The derivative of (I3.13P 
is 

f' {z) ^ -2N(f{l~ Nz)ip{z) (3.14) 

where 

1 1 - N-'^z 

^W-^p+T;^- (3.15) 

Four cases must be considered 

1° 0<iV<l, £<N-'^-l (3.16) 
2° < iV < 1, > _ 1 (3.17) 



3° iV>l, d2>i^j(8Ar2_6)yi-^^+8iV2-9j (3.18) 
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4° N>1, d^<^^(^{8N^-6)^l~^^+8N^-9j (3.19) 

Each case is illustrated in Figures 1-4. For greater convenience, the roots 
of p.l2p are shown as the intersection points of the two graphs rj = F{z) and 
7? = -A. 

The results are the following: 

In case 2 the dispersion relation defined by (|3.12p - (|3.13p has four real roots 
for A not large. 

In case 3 the dispersion relation has four real roots for " intermediate" values 
of A. 

In other cases the dispersion relation has two real roots and two complex 
roots. 

Precise formulations are given in the following Propositions. 

Proposition 1 // inequalities 113.16]} are satisfied, equation has two real 

roots and two complex roots for any positive A (Figure 1). 

Proposition 2 If inequalities {3.11^ are satisfied, the equation ip{z) — defined 
by \3.15\) has one root zqi on interval 12- Equation US. 12]) has four real roots 
for < j4 < —F{zqi) (Figure 2a) and two real and two complex roots for 
A > -F{zqi) (Figure 2b). 

Proposition 3 If inequalities Ii3.18\) are satisfied, the equation ip(z) = defined 
by i3.15]) has two roots zo2 and zqs on interval I2 ■ Equation 113. 12\) has four real 
roots for -~F{zo2) < A < — F(zo3) (Figure 3b), two real and two complex roots 
for A < -F{zo2) (Figure 3a) and A > -^(zos) (Figure 3c). 

Proposition 4 // inequalities i3.19\} are satisfied, equation i3.12]) has two real 
roots and two complex roots for any positive A (Figure 4)- 

The proofs are given in Appendix. 

It is convenient to rewrite conditions (|3.16p - (|3.19p in terms of b, k and M : 



1° : < Af < 1, fc < ^— — — (3.20) 

bM ^ ' 



r: 0<M<1, k<±^ orM>l,fc<^ (3.21) 

< fc (3.23) 



Mo(Af) 

Here do(Af) is the root of the equation f{d,M) — 0, where 
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(0 < do{M) < M). This equation has a unique solution d = do{M) in a domain 
< d < M because of the properties fd{d, M) > 0, f{d, M) — > — oo as d ^ 
and f{M,M) = AP > 0. 

InequaUties (|3.20p - (|3.23p are illustrated in Figure 5. For any given flow 
parameters M and b, one can find a sufliciently large wave number k satisfying 
inequality (|3.23p . It means that the flow is always unstable with respect to 
perturbations with sufficiently short wave lengths (as it follows from Proposition 
4, the dispersion relation has complex roots). The fiow with subsonic relative 
velocity is stabilized in an intermediate interval of wave lengths if the pure 
liquid layer is thin enough with respect to the bubbly layer. In the case of 
supersonic relative velocity a similar stabilization is observed for long waves 
(see Proposition 2 and inequalities (|3.2ip ). For another intermediate interval 
of wave lengths (|3.22[) the stabilization of perturbations is attained for layer 
depths of the same order (see Proposition 3). 

4 Conclusion and discussion 

We have derived from Hamilton's principle of stationary action governing equa- 
tions and boundary conditions at the contact interfaces in bubbly fluids. It 
has been shown that the dynamic condition on the interface reduces to the 
continuity of the average pressure. By using the boundary conditions derived, 
we have studied the Kelvin-Helmholtz instability of two superposed layers of 
a pure incompressible fluid and a bubbly fluid. We have shown that in con- 
trast to the case of two incompressible fluids when the instability develops for 
any length of perturbations (if the gravity, capillarity or compressibility are not 
taken into account), the presence of bubbles can stabilize the flow in some range 
of perturbation wave lengths. 

The stabilizing effect is due to the following reason. The development of 
classical Kelvin-Helmholtz instability leads to the appearance of wave-like bulges 
at the interface between two fluids. The presence of bubbles in a fluid (i.e. new 
interfaces) permits one to transform a part of the energy responsible for the 
bulge formation into the energy of radial oscillations of bubbles. 

The viscosity effect on the flow stabilization is an important issue. The 
method of viscous potential flows developed by Joseph et al (1999) in the analysis 
of Ray leigh- Taylor and applied by Funada & Joseph (2001) in the analysis of 
Kelvin-Helmholtz, could be used here. 
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6 Appendix 



Real roots of equation (|3.12p belong to the intervals 

h : -(1 + ^2)1/2 < 2 < /2 : 1< z < (1 + ^2)1/2 

Assuming that U20 > uiq and the wave number k is positive, we see that 
the relative Mach number M and the constant d are positive. Equation (|3.12p 
has always at least two real roots for any positive values of d, A and TV because 
on the real axis F{z) tends to —00 as \z\ tends to 1 + 0, and F{z) = for 
\z\ = (1 + d'^y^^. Additional two real roots can appear in subsets of real axis 
where F (z) changes sign. 

Notice, that the interval Ii contains exactly one root of (|3.12p because, 
obviously, F (z) < for z < — 1 (see p.l4p '). We prove below that, depending 
on the flow parameters, equation (j3.12p has one or three real roots on interval 
h. 

First, let us show that ip' (z) > on /2 if < iV < 1. If V3/2 < N < 1, then 
the derivative 

/, , 3z2_47Vz + l 
^ N{z^ 1)3 

vanishes at points zi < Z2, where 




.2,i = ^|1±a/1-^ • (6-1) 



The inequality 93 (z) > is fulfilled because in this case the roots Zi satisfy the 
inequalities < zi < Z2 < Ij and the interval where (z) < has no common 
points with I2. For < iV < \/3/2 the roots Zi of the equation (/? (z) = are 
complex and (z) > on /2. 

Proof of Proposition 1. It follows from inequalities (|3.16[) that 1 — iVz > 
for z e /2. Let us show that F (z) > on this interval. It was proved above 
that iy9 (z) > on /2 if < iV < 1. Taking into account that 



vWdP + 1) = ^^U^ (iVVrf2 + 1 - 1) < (6.2) 



we obtain inequalities <p{z) < 0, F (z) > on interval (1, v^P+T)- Then the 
equation (|3.12p has only one root on l2- Proposition 1 is proved. 

Proof of Proposition 2. In this case N^^ E I2, F {z) > for 1 < 

z < N^^ and F(z) has a local maximum at z = N^^, where F {N^^) = 0, 
F{N^^) = . Taking into account that (^'(z) > , ip{V(fTT) > and 
((p + i)iv2 - 1 

(p{N ) = — (p^pf2 — thai there exists the unique point zoi on 

I2, zoi G (A^"\ V(PTT), such that (^(zqi) = 0, i^'(zoi) = 0. At this point F(z) 
takes a local negative minimum: i^(zoi) < 0. It is obvious that equation (|3.12p 
has three roots on I2 only if < A < — i^(zoi) (see Figure 2a). Proposition 2 is 
proved. 
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Proof of Proposition 3. When inequalities (|3.18p are satisfied, the func- 
tion (p{z) is positive in a neighbourhood of the end points of I2 {^(z) — > 00 as 
z ^ 1 + 0, (p{y/d'^ + 1) > 0) and 1 - A^z < . The function F{z) can be non 
monotone only in the case when ip{z) takes negative values and, consequently, 
(f (z) vanishes at some point of l2- It follows from ()6.ip that the roots of the 
equation (p {z) = are such that zi < 1, Z2 > 1 for iV > 1. It is easy to verify 
that a stronger inequality Z2 > is valid for > 1. Hence only the root Z2 of 
the equation ip (z) = can belong to 12- The inequalities 



d2>z|-l, 6) ^(Z2) - ^ + i^^i^^ < (6.3) 



provide the inclusion Z2 G I2 and the existence of the roots Z02, Z03 of the 
equations ip{zQ2) = F (Z02) = 0, (fizos) — F (zqs) — satisfying inequalities 
1 < Z02 < zo3 < ^/<P + 1. The function F{z) has a local maximum at Z02 and 
a local minimum at zo2. It means that equation (|3.12p has three roots on I2 
if and only if — F(zo3) < A < — F(zo2) (see Figure 3a). For A > —F{zq2) and 
< A < —F(zq^) this equation has only one root (see Figures 3b, 3c). Notice 
that inequality (|6.3b ) is a consequence of (|6.3b l because for iV > 1 

>zl-l (6.4) 



N-^Z2 - 1 

Using the identity 

N-^Z2-l 27 



- 6h/ 1 ^ + 87v^ - 9 



one can show that inequality ()6.3b ) is equivalent to the second inequality in 
(j3.18|) . Proposition 3 is proved. 

Proof of Proposition 4. In accordance with (|6.4|) two cases are possible: 



Inequalities (|6.5b ,) are equivalent to the inclusion Z2 G I2 (see ()6.3p ) and the 
positiveness of minimal value of ip{z) on 12- Using the inequalities (p{z) > 0, 
1 — A^z < and (j3.14p we show that F (z) > 0. Consequently, equation (|3.12p 
has only one root on l2- 

If inequality (|6.5b ) is fulfilled, then Z2 ^ I2 and (p (z) does not change sign 
on 12 - The function (p{z) is positive on I2, and F{z) is a monotone function. We 
see that equation ()3.12|) also has only one root on l2- Proposition 4 is proved. 
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Figure 1: The variation of F as a function of z for case 1. The dispersion 
relation F{z) + A = has only two real roots z*, i = 1,2 for any positive A. 
Graph is drawn for N — 0.25 and d = 2^/2. 
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Figure 2: The variation of F as a function of z for case 2. The dispersion relation 
Fix) + A = has four real roots z*, i = 1, 2, 3, 4 for not large values of A and 
two real roots z*, i = 1, 2 for large values of A. Graph is drawn for = 3 and 
d= 7150. 
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Figure 3: The variation of _F as a function of z for case 3. The dispersion 
relation F(z) + A = has two real roots z*, i = 1, 2 for A small or large, and 
four real roots z*, i = 1, 2, 3, 4 for A intermediate. Graph is drawn for N = 1.1 
and d = 4.45. 
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Figure 4: The variation of F as a function of z for case 4. The dispersion 
relation F{z) + A = has only two real roots z*, i = 1,2 for any positive A. 
Graph is drawn for = 2 and d— 10.5 
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Figure 5: Cases 1-4 determined by inequahties (|3.20p - (|3.23|) define four con- 
nected sets in (M, c?)-plane. In regions 2~ and 3~ four real roots can exist. 
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